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Abstract. The most intriguing question related to the chemical evolution of protoplanetary 
disks is the genesis of pre-biotic organic molecules in the planet-forming zone. In this contribution 
we briefly review current observational knowledge of physical structure and chemical composition 
of disks and discuss whether organic molecules can be present in large amounts at the verge of 
planet formation. We predict that some molecules, including CO-bearing species such as Eh CO, 
can be underabundant in inner regions of accreting protoplanetary disks around low-mass stars 
due to the high-energy stellar radiation and chemical processing on dust grain surfaces. These 
theoretical predictions are further compared with high-resolution observational data and the 
limitations of current models are discussed. 
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z£ • 1. Introduction 

The origin and evolution of life as we know it are tightly related to the chemistry 
of complex carbon-bearing molecules. While the transition from macromolecules to the 
simplest living organisms has likely proceeded on Earth, we do not know yet organic 
y—{ ' molecules of what complexity have been available during build-up phase of the primor- 
dial/secondary Earth atmosphere and oceans. During the last few decades a multitude of 
species, including alcohols (e.g. CH 3 OH), ethers (e.g. CH 3 OCH 3 ), acids (e.g. HCOOH) 
have been discovered in interstellar space, with atomic masses up to a few hundrecjj 
(for a recent review see |Snyder 2006[ ). A precursor of amino acids, amino acetonitrilc, 
and the simplest sugar, glycolaldehyde, have been found toward the star-forming region 
Sagittarius B2(N) (jHollis et al. 2004j iBefloche et al. 20081) . Thus > man Y simple "blocks" 
of prebiotic molecules do exist in space. It is natural to ask what happens to these species 
during the prestellar/hot core phase, passage of an accretion shock, and inside a proto- 
planetary disk. Are organic molecules present in large amounts in circumstellar disks at 
the verge of planet formation? Could they form and survive in such a harsh environment 
as an accretion disk? 

Despite the variety of "interstellar" molecules, only formaldehyde (H2CO) and a few 
other non-organic species have been detected and spatially resolved with interferome- 
ters in several nearby protoplanetary disks (e.g., |Dutrey et al. 1997[IKastner et al. 19971 



lAikawa et al. 2003, Qi et al. 2003, Dutrey et al. 2007b). These multi-molecule, multi-transition 



studies allowed to constrain basic disk parameters like radii, masses, kinematics, tempera- 
ture and density profiles, ionization degree and depletion factors (e.g., Dartois et al. 2003, Se menov et al. 2005l Qi 

] http://astrochemistry.net/ 

1 



Th. Henning and D. Semenov 



et al. 2006,Pietu et al. 2007 Qi et al. 2008). Using advanced chemical models and indirect 



observational evidence, one can get a clue about the presence of other, yet undetected 
organic molecules in disks and estimate their abundances. 

It is commonly believed now that the formation of complex (organic) molecules starts 
already in cold dense cloud cores on dust grain surfaces serving as catalysts for many 
exothermic reactions between radicals and light atoms, with formaldehyde being one of 
the precursors for complex organic molecules. The newly produced species can be eventu- 
ally returned into the gas phase either during the slow heat-up phase after the formation 
of a central star (Garrod & Herbst 2006) or due to some thermal/non-thermal desorption 
mechanisms, like cosmic rays/X-rays/UV heating of grains (d'Hcnd ecourt et al. 1 982. 
|Leger et al. 1985| |ShaIabiea fc Greenberg 1994[|Najita et al. 2001[ lGa rrod et al. 2007) . 

Transformation of a cloud into an actively accreting disk caused by gravitational col- 
lapse and angular momentum transport further modifies the composition of gas and 
dust as it passes through a shock front (jLada 1985} IHas scl 2004). Furthermore, dy- 
namical transport in accretion disks can also be efficient for enriching the gas with 
complex species through evaporation of icy mantles in warmer, less opaque regions 
( |Willacy et al. 2006} ISemenov et al. 2 006 ) . The UV radiation from the star and inter- 
stellar radiation field plays a major role in disk chemistry by dissociating and ionizing 
molecules, and heating the gas above the midplane where many molecular lines get ex- 
cited (e.g.. lvair Zadelho ff et al. 20 03). Disk chemistry by itself can lead to the production 
of complex organic molecules. 

The recent detection of Nell line emission from several protoplanetary disks by Spitzer 



(|Pascucci et al. 20"07| supports theoretical predictions of Glassgold et al. (2007) that the 
upper disk parts can be ionized and heated by the intense X-ray radiation from a 
young star. This thermal bremsstrahlung radiation is likely produced in reconnection 
loops in stellar corona at a distance of up to 0.1 AU from the star and is capable 
of penetrating deeply into the disk inner region - the zone where planets form (e.g., 
Igea & Glassgold 1999[ ). This X-ray radiation ionizes helium atoms, which destroy CO 
and replenish thedisk gas with ionized atomic carbon. This leads to the formation of 
heavy cyanopolyynes and long carbon chains, partly on grain surfaces, which may lock 
a significant fraction of elemental carbon in the inner disk region (Semenov et al. 2004). 
Thus it is of utter importance to reveal which mechanisms and processes are impor- 
tant during various stages of protoplanetary disk evolution, using advanced theoretical 
modeling and high-quality observational data. 

In the framework of the "Chemistry In Disks" (CID) collaboration between groups in 
Heidelberg, Bordeaux, Paris and Jena, we have initiated a program to study and char- 
acterize chemical evolution and physical properties of nearby protoplanetary disks sur- 
rounding young stars of various masses and ages (see, e.g., |Dutrey et al. 2007b[ ). Among 
the various goals of the project, we searched for emission lines of precursors for complex 
organic molecules and detected and resolved the disk around the low-mass star DM Tau 
in the H2CO(3-2) line with the Plateau de Bure interferometer (though with modest 
SNR of ~ 3 — 5). We found that the H2CO emission is not centrally peaked as is the 
dust continuum, but shows an asymmetric, ring-like structure with a large inner "hole" 
of ~ 100 AU. The high-resolution observations of DM Tau with the IRAM interferome- 



ter by Dartois et al. (2003) | and Pietu et al. (2007) did not reveal the presence of central 



depression either in CO lines or in dust continuum. In this contribution we review our 
current knowledge about chemistry in disks and predict that the H2CO inner hole, if 
it really exists, should likely be caused by chemical effects, which we will discuss in the 
paper. 
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2. Observational facts 

Up to now more than 150 molecules have been discovered in spac^j]. Among them only 
a small fraction have been detected in planet-forming disks with the aid of millimeter 
interferometry: CO and its i sotopes, CN, HCN , DCN, HNC, H 2 CO, C 2 H, CS, HCO+, 
H 13 CO+, DCO+, and N 2 H+ (|Dutrey et al. 1997[lKastner et al. 1997HAikawa et al. 20031 



|Dutrey et al. 2007a Qi et al. 2008 ). A few bright and large disks like surrounding DM 
Tau, LkCa 15, and MWC 480 have been investigated in detail in a dozen of molecular 
transition and dust continuum. Combined analysis of these line and continuum data al- 
lowed us to derive disk sizes, orientation, kinematics, distribution of temperature, surface 
density, and molecular column densities (e.g., Dartois et al. 2003, Qi et al. 2006,Isc lla et al. 20 07,Pietu 
et al. 2007). 

The lines of the abundant CO molecule serve as a probe of disk geometry as well 
as thermal structure and surface density distribution and kinematics. Due to selective 
photodissociation, disks appear increasingly larger in the dust continuum and the C 18 0, 



13 CO, and 12 CO lines, with a typical radius of - 300-1000 AU (IDutrey et al. 2007aA 



An important finding is the presence of vertical temperature gradients in many disks, 
as predicted by physical models, while a few disks with large inner cavities do not show 
evidence for such a gradient (e.g., GM Aur). Furthermore, disks around hotter Herbig 
Ae/Be stars are systematically warmer than those around cool Sun-like T Tau stars. 



Recently, Qi et al. (2006) have found that the observed intensity ratios of the CO low- 
to high-level lines in the TW Hya disks require an additional heating source, which could 
be the X-ray stellar radiation. 

The observed lines of C 2 H, CN, and HCN are sensitive to the properties of the im- 
pinging UV radiation, in particular to the fraction of the total UV luminosity emit- 
ted in the Ly a line ( |Bergin et al. 200"3| ). The observed ratio of the CN to HCN column 
densities in disks is typical of photon-dominated chemistry, as predicted by the chem- 
ical models (Chapillon et al. 2008, in prep.). Molecular ions (HCO + and N 2 H + ) are 
the dominant charge carriers at intermediate disk heights and their observations al- 
lowed to constrain the ionization fraction in these regions, with a typical value of 10 -8 
QQi et al. 20031 |Dutrey et al. 2007bp . The observed ratios of DCO+ to HCO+ and DCN 
to HCN column densities have much higher D/H value than the cosmic abundance of 
~ 0.01% and thus deuterium fractionation is effective in disks ( |Qi et al. 2008| ). 

In general, the observed molecular abundances are lower by factors 5-100 in protoplan- 
etary disks compared to the values in Taurus molecular cloud, likely due to efficient freeze- 
out and photodissociation. A puzzling observational fact is that a significant reservoir of 
very cold CO and HCO + gas exists in the disks of DM Tau and LkCa 15, at temperatures 
< 13-17 K, which are below the freeze-out temperature of CO (20 K). Conventional chem- 
ical models cannot explain this fact without invoking a non-thermal desorption mech- 
anism that works in the dark disk midplane, like efficient turbulent diffusion and UV - 



photodesorption driven by cosmic ray particles (jSemenov et al. 2006|, Oberg et al. 2007 ). 

Observations of dust thermal emission at (sub-) millimeter and centimeter wavelengths 
are used to measure the slope of the wavelength dependence of the dust opacities that is a 
sensitive indicator of grain growth and sedimentation in disks (e.g., Rodmann ct al. 2006). 
There is strong evidence that in many evolved disks, with ages of a few Myrs, dust grains 
grow until at least pebble-like sizes. The results from the Infrared Space Observatory and 
Spitzer reveal the presence of a significant amount of frozen material and a rich variety of 
amorphous and crystalline silicates and PAHs in disks fe.g. Jvan den Ancker et al. 2 000, 
Ivan Dishocck 2004, B ouwman et al. 2008|) . The PAH emission features at near- and mid- 
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Figure 1. Physical and chemical structure of a protoplanetary disk. In the dark, dense and 
cold midplane most of molecules reside on dust grains, and chemical evolution is dominated 
by ion-molecule and surface reactions. This region has the lowest ionization degree, with dust 
grains being the most abundant charged species. A warmer intermediate layer is located above 
midplane. It is heated by mild UV radiation. Many reactions with barriers can occur and a rich 
variety of molecules exist in the gas phase. This is the zone where most of the molecular lines 
are formed. The ionization fraction in the intermediate layer is determined by a multitude of 
molecular ions, in particular HCO + . Further above, a hot, highly ionized disk atmosphere is 
located, where only the simplest radicals and ions (apart from H2) survive. This is the region 
where ionized carbon is abundant and C + emission lines are excited. 



infrared wavelengths are excited by the incident stellar radiation field and as such depend 
on disk vertical structure and turbulent state (jDullemond et al~2 007). These lines are 
more easily observed in disks around hot, intermediate-mass Herbig Ae/Be stars as com- 
pared to cool, Sun-like T Tauri stars (e.g. JAcke fc van den Ancker 2004|,lGeers et al. 2 007, 
|Sicilia-Aguilar et al. 2007[ ). 

Various solid-state bands observed at 10-30/mi in emission belong to amorphous and 
crystalline silicates at T > 100 — 300 K with varying Fe/Mg ratios and grain topol- 
ogy/sizes (e.g., [van Boekel et al. 20041lNatta et al. 2007llBouwman et al. 2008||Voshchinnikov fc Henning 20080 . 
The composition of the hot gas in the inner disk as traced by ro-vibrational emission lines 
from CO, CO2, C2H2, HCN and recently H 2 and OH, suggests that complex chemistry 
driven by endothermic reactions is at work there (jBrittain et al. 2003|,ILahuis et al. 2006} 
Eis ner 20*07]. |Salyk et al. 20 08 ) . At larger distances from the star the disk becomes colder 
and most of these molecules stick to dust grain, forming icy mantles. The main man- 
tle component is water ice with trace amount of other more volatile materials like CO, 
C0 2 , NH 3 , CH 4 , H 2 CO, HCOOH and CH 3 OH (jZasowski et al. 20070 . Typical relative 
abundances of these minor constituents are about 0.5-10% of that of water. 



3. Chemical structure of a disk 

The chemical evolution of protoplanetary disks has been investigated in detail by using 
robust chemical models ( |Willacy et al. 1998jlAikawa fc Herbst 19 99, Will acy fc Langer 2000] 



lAikawa et al~2 002, Bcrgin ct al. 2003, van Zadelh offet al. 2003i|llgner et al. 2004] Semenov 
et al. 2004 |Willacy et al. 20061 |Hgner fc Nelson 20080 . The current theoretical picture 
based on a steady-state prescription of the disk structure divides the disk into three 
zones, see Fig [TJ Before planets have formed and disk gas is dispersed, the dense mid- 
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plane is well shielded from stellar and interstellar high-energy radiation. While its inner 
part can be heated up by accretion, the outer zone is cold, T ~ 10 — 20 K. The only 
ionization sources are cosmic ray particles and decay of short-living radionuclides, and 
thus matter remains almost neutral, with a low degree of turbulence. The molecular com- 
plexity in the midplane is determined by ion-molecule and surface reactions, with most 
molecules sitting on the grains. Adjacent to the midplane a warmer zone is located, which 
is partly shielded from stellar and interstellar UV/X-ray radiation. The complex cycling 
between efficiently formed gas-phase molecules, accretion onto dust surfaces, rapid sur- 
face reactions, and non-negligible desorption result in a rich chemistry (see for review 
Bergin et al. 2007D . The inner part (< 10 — 20 AU) of this region in the disks around 
T Tauri stars can be substantially ionized by stellar X-ray radiation. The intermediate 
molecular layer is sufficiently dense (~ 10 5 — 10 6 cm" 3 ) to excite most of the observed 
emission lines. Atop a hot and heavily irradiated surface layer is located, where C + , light 
hydrocarbons, their ions, and other radicals like C2H and CN are able to survive. This 
is the region where PAH and silicate emission features are produced. 



4. Theoretical constrains 

Here we discuss a puzzling observation of a putative H2CO inner cavity in the disk 
of DM Tau. The Plateau de Bure interferometric image of the DM Tau disk at the 
1.5" resolution in the H2CO (3i2-2i2) line is shown in Fig. [2] (left panel). Despite high 
noise level, the H2CO emission appears as asymmetric ring-like structure, with a dip in 
southern direction. To make a proper analysis of these data, a consistent combination 
of disk physical and chemical models along with radiative transfer in molecular lines is 
used. To simulate the disk physical structure we utilize a 1+1D flared disk model which 



is similar to the model of D'Alessio et al. (1999) with a vertical temperature gradient. 



The dust grains are modeled as compact amorphous silicate spheres of uniform 0.1 /jm 
radius, with the opacity data taken from Semenov et al. (2003)| and a dust-to-gas mass 



ratio of 100. The accretion rate is assumed to be 2 x 10 y M Q yr 1 , a = 0.01, and the disk 
outer radius is 8000 AU. We focus on the observable disk structure beyond the radius of 
~ 10 AU. The total disk mass is 0.07 Mq and the disk age is 5 Myr (Pietu et al. 2007). 

We assumed that the disk is illuminated by UV radiation from the central star with 
an intensity x — 410 xo a t 100 AU and by interstellar UV radiation with intensity Xo 
in plane-parallel geometry (jDraine 19781 Ivan Dishoeck et al. 20061 Dutrey et al. 2007b[). 



We model the attenuation of cosmic rays (CRP) by Eq. (3) from Semenov et al. (2004) 
with an initial value of the ionization rate Ccrp = 1-3 • 10~ 17 s _i . In the disk interior 
ionization due to the decay of short-living radionuclides is taken into account, assuming 
an ionization rate of 6.5 • 10~ 19 s _1 (Finocchi & Gail 1997]). The X-ray ionization rate 



in a given disk region is computed according to the results of Glassgold et al. (1997a) 



Glassgold et al. (1997b) with parameters for their high-metal depletion case and a total 



X-ray luminosity of ~ 10 30 erg cm 2 s 1 ([Glassgold et al. 2005). The gas-phase reaction 



rates are taken from the RATE 06 database (jWoodall et al. 2007| , while surface reactions 



together with desorption energies were adopted from the model of Garrod & Herbst (2006) 
A standard rate approach to the surface chemistry modeling, but without H and H2 tun- 
neling was utilized (jKatz et al. 1999|) . 

Using the time-dependent chemical code "ALCHEMIC'IJ, we simulated 5 Myr years of 
evolution in the DM Tau disk, followed by 2D non-LTE line radiative transfer modeling 



with "URANIA" ( |Pavlyuchenkov et al. 20"07| , see Fig. GU The observed ring of H 2 CO 
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Figure 2. (Left) The observed integrated intensity map of H2CO (3-2) in the disk of DM 
Tau shows an asymmetric shell-like distribution with a chemical inner "hole" of ~ 100 AU 
in radius. (Right) The same features are present in the synthetic integrated intensity map of 
para-EhCO (4-3) that is produced with a realistic disk physical and chemical model and line 
radiative transfer. 

emission with a depression is fully reproduced by our model. We found two possible 
explanations why such a large-scale chemical hole can exist in the disk around DM Tau. 
The ~ 100 AU hole in the H2CO emission is fully reproduced by both a disk model 
with X-ray driven chemical processes and somewhat less markedly in the model without 
surface chemistry. These two models predict different spatial distributions of molecular 
species, which can be tested by future interferometric observations. 

Destruction of formaldehyde has important consequences for organic chemistry. The 
X-ray chemical model leads to the clearing of an inner hole of ~ 100 AU radius in all 
chemically related CO-bearing species, including HCO + , by converting gas-phase CO into 
heavier C02-containing and chain- like hydrocarbon molecules. In contrast to CO, these 
heavier species are locked on dust surfaces in the inner disk region, where temperatures 
are lower than about 35-50 K (Fig. [31 solid line). This implies substantially different 
initial conditions with respect to the presence of complex organic molecules inside the 
planet-forming zone of protoplanetary disks, if this X-ray driven chemistry is important. 

The less realistic model without surface chemistry shows the inner depression in col- 
umn densities of highly saturated molecules only, like H2O, NH3, and to some extent 
H2CO (see Fig. 02 dashed line). These species are formed on dust surfaces in a sequence 
of hydrogen addition reactions. Though at current stage we cannot fully distinguish be- 
tween these two scenarios, for our understanding of the evolution of organic species in 
protoplanetary disks it will be of great importance to verify which of these explanations 
are valid. 



5. Summary 

We briefly overview recent progress in our understanding of chemical evolution in 
protoplanetary disks, from both the theoretical and observational perspective. A puzzling 
observation of the chemical inner hole visible in the spatial distribution of the H2CO 
emission in the disk of DM Tau is addressed theoretically. We found that such a hole can 
be explained either by the absence of efficient hydrogenation reaction on dust surfaces 
or efficient processing of disk matter by stellar X-ray radiation in the inner disk region, 
which was overlooked in previous studies. In future, when the Atacama Large Millimeter 
Array will become operational, the planet-forming zone of disks will be observable and 



Organic matter in disks 



7 




r, AU r, AU 

Figure 3. (Left) The radial distribution of the HCO + column density in the DM Tau disk as 
computed with 3 different chemical models: 1) the stellar X-ray luminosity is assumed to be 
close to the observed value of 10 30 erg/s (solid line), 2) no X-ray radiation penetrates into the 
inner disk (dotted line), and 3) the model without surface reactions but with Lx = 10 30 erg/s 
(dashed line). (Right) The same calculations but for the chemically related H2CO molecule. 

this hypothesis can be verified. In general, chemo-dynamical models of disks together with 
intcrferomctric observations well lead to a comprehensive understanding of the molecular 
inventory of protoplanetary disks. 
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